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Abstract. Advances in recent years have made it possible to explore quantum dots as a viable technology for scalable 
quantum information processing. Charge qubits for example can be realized in the lowest bound states of coupled quantum 
dots and the precision control of the confinement potential allows for the realization of a full set of universal qubit 
gates, including arbitrary single-qubit rotations and two-qubit C-NOT gates. In this work we describe a novel scheme for 
implementing quantum random walks on arbitrarily complex graphs by extending these elementary operations to the control of 
a two-dimensional quantum dot grid. As single-qubit rotations constitute the essential building blocks of our implementation 
scheme, we also present numerical simulations of one such mechanism by directly solving the corresponding time-dependent 
Schrodinger equation. 
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INTRODUCTION 

The remarkable speed-up observed in the Shor factor- 
ization algorithm [1] and the Grover search algorithm 
[2] has prompt intense interest in developing appropriate 
physical systems for quantum computation and quantum 
information processing. Among the various proposals, 
solid-state systems are particularly attractive, since they 
are potentially ready to be integrated into large quantum 
networks and are enticing to the present semiconductor 
industry. 

Quantum Dots [3] have in particular inspired many 
solid-state based proposals for quantum information pro- 
cessing and are routinely manufactured in numerous lab- 
oratories. When confining exactly one electron, pairs 
of coupled quantum dots can naturally form qubits, the 
building blocks for quantum circuits. Charge qubits for 
example can be formally defined as a\L) + j3|7?) where 
a and j3 represent the single electron amplitude to be 
present in the the lowest bound states of the left and right 
quantum dot respectively. There have been a number of 
proposals for performing quantum operations on solid 
state charge qubits [4, 5, 6, 7]. In particular references 
[8, 9] have demonstrated that a precision control of the 
confinement barrier potential allows for the realization 
of arbitrary single-qubit rotations. 

In this paper we utilize these elementary qubit opera- 
tions as the basis of a novel scheme for physically imple- 
menting quantum random walks on any arbitrary undi- 
rected graph. To do so we will first provide a theoretical 
platform for quantum random walks on graphs followed 
by a description of our proposed physical implementa- 
tion scheme. Finally we provide the results of our numer- 



ical simulations for performing robust qubit rotations, 
which are fundamental for implementing of the walk. 



QUANTUM WALK ON A GRAPH 
Motivation 

Random walks have been employed in virtually every 
science related discipline to model everyday phenom- 
ena such as the DNA synapsis [10], animals' foraging 
strategies [11], diffusion and mobility in materials [12] 
and exchange rate forecast [13]. They have also found 
algorithmic applications, for example, in solving differ- 
ential equations [14], quantum monte carlo for solving 
the many body Schrodinger equation [15], optimization 
[16], clustering and classification [17], fractal theory [18] 
or even estimating the relative sizes of Google, MSN and 
Yahoo search engines [19]. 

Whilst the so called classical random walks have been 
successfully utilized in such a diverse range of applica- 
tions, quantum random walks are expected to provide us 
with a new paradigm for solving many practical prob- 
lems more efficiently [20, 21]. In fact quantum walks 
have already inspired efficient algorithms with applica- 
tions in connectivity and graph theory [22, 23], as well 
as quantum search and element distinctness [24, 25], due 
to their non-intuitive and markedly different properties, 
including faster mixing and hitting times. 

To illustrate, take Alice, who is a classical random 
walker originally positioned at x = 0. To decide whether 
to take a step to the left or right, he flips a coin with the 
two possible outcomes labeled by + and — and respec- 



tive probabilities and P-. Looking at the outcome, 
she then knows with certainty which way to move; either 
to the left {x = —I) or io the right (x = +1). Hence af- 
ter many iterations of this process Alice's classical walk 
traces a single path within a decision tree and the prob- 
ability for finding him at a given position x follows a 
Gaussian distribution. Meanwhile Bob, who is a quan- 
tum walker, flips his coin but never looks at the outcome. 
Instead he steps simultaneously to the left and right with 
complex amplitudes s^j^ and s^- such that = / + 

and \s^-\^ = P-. After many interactions Bob's quantum 
walk results in a probability wavefunction with a finite 
amplitude to be present everywhere in the tree. At the 
end we can get Bob back in one piece by "collapsing" 
his wavefunction. This will allow him to emerge at some 
position X with a probability given by the peculiar distri- 
bution depicted in Fig 1 . 
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FIGURE 1. A comparison of quantum vs. classical random 
walk. A classical random walk leads to a Gaussian distribution 
(dotted) while its quantum counterpart produces an un-intuitive 
distribution (solid) which expands very rapidly. 

Over the last few years there have been several pro- 
posals for such a physical implementation using a va- 
riety of physical systems including NMR [26], cavity 
QED [27, 28], ion traps [29], classical and quantum op- 
tics [21, 30, 31, 32, 33], optical lattice and microtraps 



[34, 35] as well as quantum dots [36, 37]. None of the 
existing proposals however consider quantum random 
walks on general graphs, with the majority describing 
only a one-dimensional implementation. This is while 
from an application point of view most useful algorithms 
would involve traversing graphs with arbitrarily complex 
structures. In this paper, we build on an earlier imple- 
mentation scheme based on optical lattices [38], and ad- 
vance a proposal for realizing a universal quantum ran- 
dom walk capable of traversing any general undirected 
graph using a two-dimensional grid of quantum dots. 



Theory 

Let us first consider a complete graph with all pos- 
sible connections between the jV nodes including self 
loops. We will relax this constraint later by removing 
the unwanted connections. Here the walker requires an 
sided coin for moving from one node to jV other 
nodes. The complete state of the walker is therefore de- 
scribed by \\if) =LjfiL^i^-,y^|j,^), where |j) denotes 
the nodes or position states of the walker, |^) specifies 
the state of the coin, and is a complex amplitude. A 
coin flip in the context of this quantum walk corresponds 
to a unitary rotation of the coin states at every node j 
using an ^ X ^ matrix Cj also known as the coin oper- 
ator. The coin operation is followed by the walker step- 
ping from node j simultaneously to all other nodes on 
the graph using a conditional translation operator f such 
that fly, ^) — > \ according to some predefined rule, 
where j and / label the two nodes at the end of an edge 
^jf [39]. The quantum walk evolves via repeated appli- 
cations of the coin followed by the translation operator. 
More explicitly 

where Ii/Aq) is the initial state of the walker, ly/n) corre- 
sponds to the state of the walk after n steps, x 
matrices % and 7) are the coin and translation operators 
at the ith step, and ^ incorporates the individual coin 
operators ci . . . which simultaneously act on all the 
nodes. The operators c can in principle invoke different 
rotations at each node j, but are often uniformly set to be 
the Hadamard matrix. 

In [38] we have shown that for any arbitrary graph 
the Hilbert space of the walk can in fact be represented 
as a two-dimensional ^ x ^ grid (Fig 2) where each 
individual node on the graph corresponds to a grid row 
(column) and the coin states within that node are the 
individual grid elements along that row (column). There, 
it is also shown that the quantum walk evolution given in 
Eq. 1 can effectively be reduced to 

\¥n)='^n'^n-l--- "^2 ^ I V^o) , (2) 



where C^^) correspond to the appHcation of the coin 
operator for the ith step, on the grid elements that are 
grouped horizontally (vertically) as depicted in Fig. 2. 

Any general graph can now be constructed from 
its corresponding complete graph by removing the un- 
wanted edges as depicted in Fig. 2. Eliminating an edge 
ejjf corresponds to eliminating the connection between 
two states and |/, j). However instead of remov- 

ing these states from the Hilbert space, the action of the 
coin operators cf (cj) can be altered in such a way as 
to isolate the unwanted states from interacting with other 
states. 
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FIGURE 2. a) A complete graph with = 6 nodes, b) The 
quantum walk Hilbert space viewed as a 2D array, where the 
element {j,k) represents the state k under node j. A generalized 
graph can be constructed by removing edges (dotted lines) 
from a complete graph. This can be physically achieved by 
modifying the coin operators in such a way as to isolate the 
unwanted states (dotted circles) from interacting with other 
states. 



PHYSICAL IMPLEMENTATION 

We propose utilizing a two-dimensional grid of quan- 
tum dots to represent the quantum walk's Hilbert space, 



where the initial distribution of a single electron wave- 
funcion among the dots represents the initial state of the 
walk. As we will see this scheme necessitates the con- 
struction of a 2yK X 2^ grid, where every second row 
(column) of dots will be used as temporary "register". In 
what follows we present a mechanism for manipulating 
the electronic wave function throughout the grid by intro- 
ducing appropriate quantum dot interactions in a manner 
which corresponds exactly to the quantum random walk 
evolution described in Eg. 2. Naturally the resulting elec- 
tron probability wave distribution gives the final state of 
the quantum walk. 

Considering the evolution of the walk in Eq. 2, what 
we require is a means by which we can implement the ac- 
tion an arbitrary ^ -level unitary operator cJ (cJ) on the 
^ nodes along the jth row (column) of the grid. In [38] 
we showed that this can be achieved via a CS decompo- 
sition [40] which essentially reduces a general ^-level 
rotation matrix to a series of pair- wise or qubit rotations 
which can in principle be readily implemented. What 
makes this implementation scheme non-trivial however 
is the fact that the resulting pair- wise interactions are not 
limited to neighboring nodes. More precisely it can be 
shown that for a general rotation of column j for exam- 
ple, we have 

c]= n ^•(^)' (3) 

for J G [2, 4, . . . ^/2] , where the action of each ^ ( J) on 
the nodes of row j consists of ^/2 simultaneous pair- 
wise interactions between nodes kd-\-r and kd-\-r-\-d/2 
where k = 0,...,^/J- 1, d = 2,4,... ^/2 and r = 
1,2- ••J/2. Clearly for all J 7^ 2, interactions are non- 
neighboring but follow a systematic form. 

In the following sections we describe two mecha- 
nisms: 1) How to implement a pair- wise interaction be- 
tween any two neighboring quantum dots, and 2) How to 
extend this to non-neighboring quantum dots. 



Pair- wise Interactions 

A pair of neighboring quantum dots can be made to 
undergo a unitary rotation R by the precision control of 
the potential barrier between them [8, 9]. Assuming that 
the confined electron is initially in the left dot for exam- 
ple (Fig. 3), by dropping the barrier the electron wave 
packet is free to move between the two dots. Returning 
the barrier to its initial state at a precise moment in time 
makes it possible to recapture the electron but this time 
with the desired distribution across the two dots. We have 
presented a numerical simulation of this scheme in the 
last section. 

Here we adopt the usual terminology to describe a 
71 rotation as one in which the electron wave-packet is 



transferred entirely from one dot to its neighboring dot. 
This is also equivalent to a SWAP or NOT gate. Likewise 
a 7r/2 rotation represents a 50 — 50 split of the wave- 
packet initially confined to either of the dots. 

^) ^k^\J' \kr\kf 
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FIGURE 3. a) Arbitrary qubit rotations can be performed by 
a precision control of the central potential barrier, b) The barrier 
raised in time to implement a n rotation, c) the barrier is raised 
in time to implement a 7r/2 rotation. 



Quantum Dot Conveyor Belt 

To interact a pair of non-neighboring quantum dots we 
are effectively forced to move the two dots close to each 
other, apply the desired rotation and then return them to 
their original location. 

What makes this process systematic and practically 
viable is the specific pattern of pairwise interactions 
arising from Eq. 3. As depicted in Fig. 4, to implement 
the action of each ^i{d) we take the following five steps: 

1. Apply ^ /2 simultaneous 71 rotations to pairs of 
quantum dots at kd -\- r and their adjacent register 
dot. This has the effect of transferring the electron 
wave packets to the register dots. 

2. Adiabatically move the register quantum dots, 
much like a "conveyor belt" carrying the electron 
wavepackets along the register row (column). 
This can be archived by carefully designing and 
manipulating the voltage applied to the electrodes 
such that the confinement potentials experience an 
effective motion. Moving the register dots by an 
amount equal to d times the quantum dot width 
would effectively allow the amplitudes dXkd^r to 
be coupled with the amplitudes at M + r + <i/2. 

3. Simultaneously apply general rotations R to 
the newly coupled quantum dot pairs. 

4. Return the register qubits to their original location 
by reversing the adiabatic motion; 

5. Introduce another n rotation to move the amplitudes 
from register dots back to their original positions. 

We emphasis that the above steps can be carried out si- 
multaneously throughout the entire grid and the rotation 
R in step 3 can be different for each quantum dot pair. 
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FIGURE 4. The initial electron wavefunction inside the 
quantum dot grid followed by 5 steps required to implement 
the operator (J = 4). 



Modeling Qubit Rotations 

One of the key issues in quantum computation is to de- 
sign a time-dependent Hamiltonian as precisely as pos- 
sible, so that one can drive the system undergo the re- 
quired logic gate manipulation. Since the intrinsic par- 
allelism of a quantum computer comes from superposi- 
tion and entanglement where phases and amplitudes play 
an eminently important role, such a design can only be 
achieved through detailed theoretical calculations includ- 
ing all perceptible interactions in the system. 

The Schrodinger equation governing the quantum dy- 
namical evolution of few-electron systems reads 



2m* 



-y(r,) 



4718 



i>j=l '^J 



The above equation can be solved efficiently using a 
Chebychev polynomial expansion [41, 8] 

exp(-/(^niax+^min>/2)^4(a)7;(-/7?)VA(r,0), 

where £^max and £^niin are the upper and lower bounds on 
the energies sampled by the numerical grid, /^(a) are 
Bessel functions of the first kind, Tn are the Chebyshev 
polynomials. The normalized Hamiltonian is defined as 

H {2H - Emsix - Emin) /(^'max - £^min) to cnsure Con- 
vergence. 

Using such a theoretical framework, all electrons in 
the system are treated on equal footing, and they evolve 
coherently in time under the influence of each other as 
well as external fields. Fig. 5 illustrates the time evo- 
lution of an electron confined in a coupled quantum 
dot system provided by the above equation. All possi- 
ble qubit rotations in the Bloch sphere can be accom- 
plished by controlling the central potential barrier and/or 
electron-electron interactions, as demonstrated in Fig. 6 
where the amplitudes of and the phase between the |L) 
and \R) states are shown as time varies. 




FIGURE 5. Time evolution of system wave-function in cou- 
pled quantum dot system. 
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FIGURE 6. Qubit rotation in the Bloch sphere as time varies. 

CONCLUSION 

In this paper, we provided a practical recipe for effi- 
ciently implementing quantum random walks on arbi- 
trarily complex graphs using electrons trapped in a 2D 
array of quantum dots. The proposed scheme is partic- 
ularly elegant since the walker is not required to physi- 
cally step between the vertices of the associated graph, 
giving rise to a significant advantage over other existing 
schemes. Also presented is a detailed simulation of con- 
trolled qubit rotation in the Bloch sphere. Our numeri- 
cal results show that even the simplest gate operations 
induce rather intricate quantum dynamics. Nonetheless, 
with precision control of the electrodes that define the 
coupled quantum dots, arbitrary qubit rotations in the 
Bloch sphere can be accomplished. 
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